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ABSTRACT 

Using a sample of 21 damped Lyman alpha systems (DLAs) and 35 sub-DLAs, we evalu- 
ate the D-index = x 1000 from m § h resolution spectra of the Mgn A 2796 profile. 
This sample represents an increase in sub-DLA D-index statistics by a factor of four over the 
sample used by Ellison (2006). We investigate various techniques to define the velocity spread 
(Av) of the Mg ii line to determine an optimal D-index for the identification of DLAs. The suc- 
cess rate of DLA identification is 50 - 55%, depending on the velocity limits used, improving 
by a few percent when the column density of Fe n is included in the D-index calculation. We 
recommend the set of parameters that is judged to be most robust, have a combination of high 
DLA identification rate (57%) and low DLA miss rate (6%) and most cleanly separate the 
DLAs and sub-DLAs (Kolmogorov-Smirnov probability 0.5%). These statistics demonstrate 
that the D-index is the most efficient technique for selecting low redshift DLA candidates: 
65% more efficient than selecting DLAs based on the equivalent widths of Mgn and Fen 
alone. We also investigate the effect of resolution on determining the N(H i) of sub-DLAs. We 
convolve echelle spectra of sub-DLA Lyo- profiles with Gaussians typical of the spectral reso- 
lution of instruments on the Hubble Space Telescope and compare the best fit N(H i) values at 
both resolutions. We find that the fitted H i column density is systematically over-estimated by 
~ 0.1 dex in the moderate resolution spectra compared to the best fits to the original echelle 
spectra. This offset is due to blending of nearby Lya clouds that are included in the damping 
wing fit at low resolution. 
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1 INTRODUCTION 

The near-UV is a critical wavelength regime for quasar absorp- 
tion line spectroscopy. Whilst blue-sensitive ground-based instru- 
ments such as UVES at the VLT have opened the door for such 
efforts as measuring the molecular content of DLAs (Ledoux et 
al. 2003; Noterdaeme et al. 2008) and the study of the low red- 
shift Lyo- forest (e.g. Kim et al. 2002), studying neutral hydrogen 
at z < 1.6 requires a space telescope. Moreover, due to the de- 
clining number density of galactic scale absorbers, such as damped 
Lyman alpha systems (DLAs), at low redshifts, blind searches for 
these galaxies place unrealistic demands on space resources. There- 
fore, whilst the number of known DLAs at z > 1.7 is now in ex- 
cess of 1000, thanks to trawling large ground-based optical sur- 
veys (http://www.ucolick.org/~xavier/SDSSDLA/index.htmll, the 
head-count of low redshift DLAs is only around 5% of the high 
z tally (e.g. Rao, Turnshek & Nestor 2006). Characterising the low- 



* Based on observations made with ESO Telescopes at the Paranal Obser- 
vatories under programmes 078.A-0003(A) and 080.A-0014(B). 



to-intermediate redshift DLA population is important not just for 
the large fraction of cosmic time that it represents, but also because 
the detection of galactic counterparts for the absorbers is largely 
only feasible at z-^bs < L 

In order to circumvent the high cost of a blind survey for low z 
DLAs, the practice in recent years has been to pre-select DLA can- 
didates based on the detection of strong metal lines in ground-based 
spectra (e.g. Rao & Turnshek 2000). Although there is no tight cor- 
relation between the Af(Hi) and Mgn equivalent width (EW), there 
is a statistical correlation for large samples, e.g. Bouche (2008) and 
Menard & Chelouche (2008). In the most recent survey of low z 
DLAs selected by metal lines, Rao et al. (2006) found that 36% of 
absorbers with rest frame EWs of Mg n A 2796 and Fe n A 2600 > 
0.5 A were DLAs. Including the additional criterion that the EW 
of Mg i A > 0. 1 A, increases success rate for identifying DLAs 
to 42%. The remaining absorbers were found to have 18.0 < log 
yV(Hi) < 20.3 (Table 1 of Rao et al. 2006), thus spanning the range 
from Lyman limit systems (LLS) to sub-DLAs. Whereas DLAs 
have neutral hydrogen column densities A'(Hi) > 2 x 10 20 cnT 2 , 
sub-DLAs are usually defined (e.g. Peroux et al. 2003a) to have 
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19.0 < log Af(Hi) < 20.3. The upper bound of this classification 
corresponds to the transition to classical DLAs, whilst the lower 
limit is driven by the clarity of the Lya damping wing necessary 
to produce a reliable fit. The LLS and sub-DLAs are usually ex- 
cluded from the statistical calculation of quantities such as f2 eas 
(the mass density of neutral gas). The contribution of sub-DLAs to 
fl gas , and even the validity of including the sub-DLAs in the cen- 
sus for neutral gas are still hotly debated topics (e.g. Peroux et al. 
2003b; Prochaska, Herbert-Fort & Wolfe 2005). Nonetheless, the 
sub-DLAs are emerging as an interesting field of research in their 
own right, and for statistical purposes it is often desirable to sep- 
arate the DLAs and sub-DLAs (e.g. Turnshek & Rao 2002). For 
these reasons, it is desirable to develop an empirical tool that al- 
lows an observer to pre-select candidate absorption systems that 
are most suitable for their purpose. Therefore, whilst the Mg n + 
Fe ii EW selection has certainly been the key to identifying the vast 
majority of low redshift DLAs, a more robust separation of DLAs 
and sub-DLAs from ground-based spectra would be welcome. 

Ellison (2006) proposed a new way to screen Mg n absorbers 
as potential DLAs, defining the D-index as the ratio between Mg n 
EW and velocity width (see Section [3~2l . For a sample of 27 ab- 
sorbers, the success rate of the D-index in selecting DLAs was 
found to be more than a factor of two improvement over the use 
of EW alone. However, only eight of the 27 absorbers were sub- 
DLAs, and yet these lower column density absorbers are more 
abundant than their high A'(Hi) cousins (e.g. Peroux et al. 2003b). 
In this work, we re-visit the assessment of the D-index as a tool for 
the pre-selection of DLAs based on high resolution spectroscopy 
of Mg ii absorbers with an enlarged sample of 56 absorbers. 



2 SAMPLE SELECTION 



We selected sub-DLAs from the compilation of Rao et al. 
(2006) for follow-up spectroscopy at high resolution, obtaining 
new echelle spectra of 13 QSOs with 17 absorbers along their 
lines of sight. Nine of these QSOs were observed with the UV 
and Visual Echelle Spectrograph (UVES) at the Very Large Tele- 
scope (VLT) and four with High Resolution Echelle Spectrograph 
(HIRES) at the Keck telescope. The observations and data re- 
duction are described in detail in Dessauges-Zavadsky et al. (in 
preparation). In Table [TJ we summarise the details of these new 
data; emission redshifts and V-band magnitudes are taken from 
Rao et al. (2006). The new data were reduced using the pub- 
lically available pipelines UVES.popler (e.g. Zych et al. 2008, 
see also http://astronomy.swin.edu.au/~mmurphy/UVES_popler) 
and xidl HIRES redux (e.g. Bernstein et al. in preparation, 
see also http://www.ucolick.org/~xavier/HIRedux/index.html I. We 
have also added to the archival sample, thanks to the donation 
of spectra that have appeared in Meiring et al. (2006, 2007) and 
Churchill et al. (1999, 2003b). Finally, one additional spectrum 
(Q02 16+080) has been obtained from the UVES archive and re- 
reduced by us (Dessauges-Zavadsky et al. in prep.). Our final sam- 
ple consists of 56 absorbers, of which 21 are DLAs and 35 are 
sub-DLAs. We have therefore doubled the total sample size of ab- 
sorbers since Ellison (2006) and increased the number of sub-DLAs 
by more than a factor of four. In Table [2] we list the full absorber 
sample, absorber redshifts, H i and Fe n column densities and refer- 



ences for these quantitiefl Out of the 56 absorbers in our sample, 
49 have available N(Fen) column densities, of which 19 are DLAs 
and 30 are sub-DLAs. 



3 RESULTS 

3.1 Sub-DLAs at moderate resolution 

Most of the sub-DLAs in our sample have been drawn from the 
surveys of Rao & Turnshek (2000) and Rao et al. (2006), where 
the Af(Hi) column densities have been determined from moder- 
ate resolution HST spectra. Although the damping wings of sub- 
DLAs should theoretically be present in data with resolutions be- 
low that of echelle spectrographs (typically FWHM ~ 6 km s" 1 , 
or 0.1 A), there are a number of practical problems which may af- 
fect the Lyo- fit at all resolutions, including blending of nearby Lya 
forest lines, the accuracy of the continuum fit and correct sky sub- 
traction. Therefore, although random errors can be quoted for the 
iV(Hi) determinations (typically 0.05 - 0.10 dex), it is important to 
also consider any systematic effects. This is particularly relevent for 
a technique that aims to distinguish DLAs and sub-DLAs, since we 
must be confident that the absorbers are correctly classified. Indeed, 
Meiring et al. (2008) have previously suggested that low resolution 
spectra might lead to over-estimates of N(H i) when only single ab- 
sorption components are fitted. Here, we quantitatively investigate 
this possibility and assess its impact on the current study. 

The HST spectra used by Rao & Turnshek (2000) and Rao et 
al. (2006) have been compiled from various surveys and archival 
data from the Faint Object Spectrograph (FOS) and Space Tele- 
scope Imaging Spectrograph (STIS; both MAMA and CCD ob- 
servations). Although the data vary in characteristics, the typical 
FWHM resolution ~ 5 A. In order to investigate the presence of 
systematic uncertainties in the N(H i) determinations of sub-DLAs, 
we have taken the UVES spectra of the 12 sub-DLAs published 
by Dessauges-Zavadsky et al. (2003) and convolved them with a 
gaussian profile of FWHM=5A to simulate the effect of HST spec- 
tral resolution. The process of convolution effectively smooths the 
noise out of the UVES data. However, we do not re-insert any noise 
characteristics, so that our fitting tests are driven by the effects of 
resolution, not noise. We use VPFIT0 to determine the H i column 
density of absorbers in the convolved spectra. In Table [3] we give 
the original UVES best fit N(H i) values from Dessauges-Zavadsky 
et al. (2003) and our VPFIT values. The comparison is shown visu- 
ally in Figure [TJ 

The main result of this test is that we find a systematic offset 
between ^V(Hi) determined from the high resolution UVES data 
and the spectra that have been convolved to mimic HST resolution. 
The low resolution fits yield N(H i) values that are typically higher 
than the UVES-determined values by 0. 1 dex, although the discrep- 
ancy can be as large as 0.3 dex. The magnitude of the discrepancy 
appears to be driven by the amount of local Lyo- absorption. At 
high resolution, it is relatively straightforward to distinguish low 
column density Lya forest clouds from the main sub-DLA absorp- 
tion. However, at low resolution, these components become indis- 
tinguishable and the fit is driven to higher values of N(H i) to ac- 
count for the increased equivalent width. 

This effect is demonstrated in Figure [2] where we show two 

1 In cases where errors on N(Fe n) are not available in the literature , we 
assign a value of 0.05 dex. 

2 http://www.ast.cam.ac.uk/~rfc/vpfit.html. 
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Table 1. New QSO observations 



QSO 


^eiii 


V 


Instrument 


Observation date 


Exp. time (s) 


S/N per pixel 


Q0009-016 


1.998 


17.6 


UVES 


Nov. 2006 


6000 


25 


Q002 1+0043 


1.245 


17.7 


UVES 


Nov. 2006 


9000 


20 


Q0157-0048 


1.548 


17.9 


UVES 


Oct 


. 2006 


6000 


20 


Q0352-275 


2.823 


17.9 


UVES 


Oct 


. & Nov. 2006 


10240 


35 


Q0424-131 


2.166 


17.5 


UVES 


Oct 


. & Nov. 2006 


6000 


30 


Q1009-0026 


1.244 


17.4 


UVES 


Jan 


. 2007 


6000 


30 


Q1028-0100 


1.531 


18.2 


UVES 


Feb 


.2007 


10240 


20 


Q 1054-0020 


1.021 


18.3 


UVES 


Jan 


. & Feb. 2007 


10240 


35 


Q 1327-206 


1.165 


17.0 


UVES 


Feb 


. 2008 


4460 


25 


Q1525+0026 


0.801 


17.0 


HIRES 


Jul. 


2006 


1200 


10 


Q2048+196 


2.367 


18.5 


HIRES 


Jul. 


2006 


3655 


20 


Q2328+0022 


1.308 


17.9 


HIRES 


Jul. 


2006 


2500 


10 


Q2352-0028 


1.628 


18.2 


HIRES 


Jul. 


2006 


5400 


20 



extreme cases, Ql 101-264 and Q2344+0342. The former of these 
sub-DLAs appears isolated from Lyor forest blending in the UVES 
spectra and has a very well defined continuum and consequently 
an excellent H i fit. Conversely, Q2344+0342 is highly blended and 
Dessauges-Zavadsky et al (2003) found it necessary to simultane- 
ously fit a second absorption system in the blue wing of the sub- 
DLA in order to obtain a satisfactory fit. In Figure[2]we show both 
the original UVES data (top panels) and convolved data (lower pan- 
els) for the two QSOs. In each panel, the blue dashed lines show 
the profiles (at the appropriate resolution for the data shown) for 
the N(Ri) values derived from the UVES data and the red dot- 
ted lines show the values determined from VPFIT to the convolved 
data. Therefore, although the N(H i) values for the top and bottom 
panels are identical for a given sub-DLA, the resolution of the pro- 
file is matched to the effective instrumental resolution of the data 
shown. For convenience, the ,/V(H i) values are given (with the ap- 
propriate colour code) at the top of each panel. Since redshift is also 
a free parameter in VPFIT, we also show its best-fit value, as well 
as the value obtained from the UVES data (taken from Dessauges- 
Zavadsky et al. 2003). In the case of Ql 101-264, it can be seen 
that the VPFIT N(H i) column density and redshift are in excellent 
agreement with the UVES values. The same is not true of the sub- 
DLA towards Q2344+0342. From the bottom-right panel, we see 
that the UVES N(H i) value (blue dashed line) is a poor fit to the 
low resolution data. The VPFIT solution (red dotted line) is driven 
to fit the blue wing, which, in reality, is severely blended with ad- 
ditional Lya absorption. This also causes the redshift of the two 
solutions to differ by 120 km s . In the top-right panel, it can be 
seen how the VPFIT low resolution fit is a poor representation of 
the high resolution data. 



Although we have demonstrated that there is a tendency to 
over-estimate N(H i) from low resolution spectra, the effect, in gen- 
eral, is relatively small, usually ~ 0.1 dex (see Table [3}. Abun- 
dances of sub-DLAs may therefore have been previously under- 
estimated by 0.1 - 0.3 dex, and there are also implications for the 
number densities of absorbers with modest HI column densities. At 
low redshift, this effect may be mediated somewhat by the lower 
line number density as the forest thins. On the other hand, it has 
been shown that sub-DLAs are often found near (in velocity space) 
to other strong absorbers (e.g. Ellison & Lopez 2001; Peroux et al. 
2003a; Ellison et al. 2007). Statistics that rely on weighted, or total 
N(H i) values should not be greatly affected by this systematic error, 



Table 3. Resolution test: N(H i) fitting results 



QSO 



UVES N(B i) FWHM=5A N(H i) 



Ql 101-264 


1.838 


19.50±0.05 


19.51+0.05 


Q1223+175 


2.466 


19.32±0.15 


19.60+0.10 


Q 1409+095 


2.668 


19.75±0.1() 


19.83+0.05 


Q1444+014 


2.087 


20.18±0.10 


20.23+0.05 


Q1451+123 


2.255 


20.30±0.15 


20.45+0.05 


Q151 1+090 


2.088 


19.47 ±0.10 


19.80+0.05 


Q2059-360 


2.507 


20.21±0.10 


20.30+0.05 


Q21 16-358 


1.996 


20.06±0.10 


20.15+0.05 


Q2155+1358 


3.142 


19.94±0.10 


20.03+0.05 




3.565 


19.37±0.15 


19.70+0.10 




4.212 


19.61+0.10 


19.58+0.05 


Q2344+0342 


3.882 


19.50±0.10 


19.73+0.08 



since the fractional contribution of blended Lyo- forest absorption 
will have only a minor impact on the high column density DLAs 
that dominate such quantities. 

Finally, we note that the systematic effect investigated here 
does not include errors in the continuum fit. The error in contin- 
uum fitting is difficult to quantify for absorbers in general, given 
that different authors use very different techniques, and that these 
techniques themselves are often adapted to the resolution and S/N 
of the data. Errors in the continuum are often estimated to add a 
further 0.1 dex uncertainty to the determination of A'(Hi). 



3.2 The D-index 

Ellison (2006) defined the D-index to be 
EW(A) 



D 



Av(km/s) 



x 1000, 



(1) 



where the EW is that of the Mg n A 2796 line and Av is the veloc- 
ity spread of the same line. The limits of Av were determined by 
excluding 'detached' absorption components where the continuum 
is recovered, and only include the complex with the largest EW. 
The limits of /l m i n . max (the minimum and maximum wavelengths 
over which the EW is calculated) were set where the absorption 
becomes significant at the 3 cr level below the continuum. 

With a larger sample, we can experiment further with the def- 
inition of D, in order to optimize the distinction between DLAs 
and sub-DLAs. We begin by experimenting with different ways of 
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Table 2. Full absorber sample for D-index calculations 



QSO 




logAf(Hi) 


N(Fe n) 


JV(H i) reference 


Fe ii reference 




00002+05 1 


0.851 


19.08 


+ 


0.04 


1 3 SQ-i-0 04 


Ran Turnshek Hi Nestor (7006^ 


Miirnnv nnnnhli thprl 

1V1LI1 1711 V , Lill 17 LI 1711 AllCU 




Q0009-016 


1.386 


20.26 


± 


0.02 


14.32±0.04 


Rao, Turnshek & Nestor (2006) 


Dessauges-Zavadsky et al 


in prep 




520 


19.54 


+ 


0.03 


1 ^ 1 7+0 04 

1 J. 1 / 


Ran Tiirrrihpk Rr l\Tp<;tnr T900M 

I\aU, ILIlllSllCri. OL 1"IC&LU1 \Z.\J\JVJJ 


Tip ccnitTP c«7^\^in ck\7 pt cil 
l_/CSSd.LlgCS Zjdvd.Llftri.y CI til. 


in nt'i'n 
111 U1CU. 


00(12 1 +004'? 


0.940 


19.38 


+ 


0.15 


1 4 62+0 04 


Ran Tnrnshplc Rr Np^tnr r7006^ 

I v iiv t . i lii ii jiiciv oc \ jLKjyjyj j 


T~)pccan opc~7n vari cLrv pt nl 

l_/CC7 0uLlgCC7 Zjd V uUAJV V rjL til. 


in nrpn 

111 Ul^lJ. 


Q0058+019 


0.613 


20.08 


± 


0.20 


15.24±0.06 


Pettini et al. (2000) 


Pettini et al. (2000) 




O01 00-4- 1 30 


2 309 


21.37 


± 


0.08 


1 <; OQ+0 01 


Mpcc 011 frpc - / ivnrl c w\f pt n 1 ( 7 004 i 
IJCaaaUgCa ZjdvdUolvy CL til. \Z.\j\j 1 t) 


Tipc c in npc-7^\fnMC pt nl 
L/CftSdUgCS ZjdVaUSIvy CL al. 


(2004) 


ooi 17+711 


0.576 


19.15 


+ 


0.07 




Ran Tnrrwhpk Rr Np^tnr I^OOM 

1 V ilV / , ± Lll HjIICjY OC l^t^LUl \£*\J\J\J j 






Q0 157-0048 


1.416 


19.90 


+ 


0.07 


14.57+0.03 


Rao, Turnshek & Nestor (2006) 


Dessauges-Zavadsky et al 


in prep 


Q02 16+080 


1.769 


20.00 


± 


0.20 


14.51+0.09 


Lu et al. (1996) 


Dessauges-Zavadsky et al 


in prep 


\/ \J J) *J jL* L. t 


1.405 


20.18 


± 


0.18 


1 5 1 0+0 03 


Rao Turnshek Hi Nestor ("7006^ 


T~)pccan opc-7fivnric.k"v pt nl 

!_/CcjCitlLlgCcj Z~i& V uLlAJV V CL dl. 


in nrpn 

111 L71CI7. 


O0424- 1 3 1 


1.408 


19.04 


± 


0.04 


1 1 45+0 07 


Ran Tiirnshpk Rr Np^tnr <?00f\) 


T^pccan oPG-7nvanckv pt nl 

!_/CcjCitlLlgCcj Z~i& V uliciJV V CL til. 


in nrpn 

111 L71CL7. 


00454-220 


0.474 


19.45 


+ 


0.03 




Ran Tnrnshpk # Np^tnr ^OOM 

IV ilV / , XLlllljll^R. OC li^jLVJl \£*\J\J\J j 






Q0512-333A 


0.931 


20.49 


± 


0.08 


14.47+0.06 


Lopez et al. (2005) 


Lopez et al. (2005) 




Q0512-333B 


0.931 


20.47 


± 


0.08 


>14.65 


Lopez et al. (2005) 


Lopez et al. (2005) 




00823-223 


0.91 1 


19.04 


± 


0.04 


i 3 57 + o 04 

1 _7 ..J / iV7 . 


Rao Turnshek Hi Nestor (2006^ 


Meiring et al. (2007) 




00877+24 


0.525 


20.30 


+ 


0.05 


14 59+0 07 

1 *"T..J y V7 . V7 i_ 


Ran Rr Tnrnshpk ("70001 


Meiring et al. (2006) 




O0841 +1 29 

Wvu^ ill jL. y 


2.375 


20.99 


± 


0.08 


14 76+0 01 


T~lpc;fiii cip';-7n\jfinQli"v pt nl ^700n ^ 

lyCAauLlcCft ijuvaLl&IVY tTL CLi « ^ Z-AJ VJ y 


Dessauges-Zavadsky et al. 


(2006) 


Q0957+561A 


1.391 


20.30 


± 


0.10 


14.46+0.05 


Churchill et al. (2003a) 


Churchill et al. (2003a) 




Q0957+561B 


1.391 


19.90 


± 


0.10 


14.34+0.05 


Churchill et al. (2003a) 


Churchill et al. (2003a) 




Q 1009-0026 


0.840 


20.20 


± 


0.06 


14.48+0.02 


Rao, Turnshek & Nestor (2006) 


Dessauges-Zavadsky et al. 


in prep. 


01009-0026 


0.880 


19.48 


+ 


0.08 


14 37+0 07 


Rao Turnshek Hi Nestor (2006^ 


Dessauges-Zavadsky et al. 


in prep. 


101 0-0047 

1 \7 1 vy UU 7 / 


1.327 


19.81 


± 


0.07 


14 53+0 03 

i*t..7->^\7.\7j 


Ran Tiirnshplc Rr Np^tnr ("7006^ 


Meiring et al. (2007) 




O1078-01 00 


0.709 


20.04 


± 


0.07 


1 5 1 0+0 03 


Ran Tiirn\hplc ^ Np^tnr ("7006^ 


Dessauges-Zavadsky et al. 


in prep. 


Q 1054-0020 


0.951 


19.28 


+ 


0.02 


13.71+0.01 


Rao, Turnshek & Nestor (2006) 


Dessauges-Zavadsky et al. 


in prep. 


Ql 101-264 


1.838 


19.50 


± 


0.05 


13.51 ±0.02 


Dessauges-Zavadsky et al. (2003) 


Dessauges-Zavadsky et al. 


(2003) 


Q1104-180A 


1.662 


20.85 


± 


0.01 


14.77±0.02 


Lopez et al (1999) 


Lopez et al (1999) 




Ol 122-168 


0.682 


20.45 


+ 


0.05 


14 55 + o 01 


T prlmiv Rproprnn Rr Ppfitipnn f70071 

l^t^LIVJLIA., IJClgClUll OC 1 CllLl^ull (i-UUi/ 


Ledoux, Bergeron & Petitjean (2002) 


Ol 1 51 +068 


1.774 


21.30 


+ 


0.08 




Dessauges-Za.vad.sky unpublished 






Q1157+014 


1.944 


21.60 


± 


0.10 


15.46 ±0.02 


Dessauges-Zavadsky et al (2006) 


Dessauges-Zavadsky et al. 


(2006) 


Ol 706+459 


0.928 


19.04 


± 


0.04 


1 2 95+0 02 


Ran Tiirnshpk Rr Np^tnr ("7006^ 


Murphy, unpublished 




01710+17^ 


1 892 


20.63 


± 


0.08 


IS 01 +0 03 


Mpcc on (Tfic: _ x n\7cifi clr\/ pt 1 1 ( 7 OOri i 
UCftadLlgCft ZjtlVtlLlSlvy CL til. ^ Z.LfVJUy 


Dessauges-Zavadsky et al. 


(2006) 


01213-002 


1.554 


19.56 


+ 


0.02 


1 4 44+0 1 


Ran Tnrnshplc Rr Np^tnr ^7006^ 
i\ cw / , xLiiii^iidv oc v ^vywvj / 


Murphy, unpublished 




Q1223+175 


2.466 


21.44 


+ 


0.08 


15.16+0.02 


Dessauges-Zavadsky unpublished 


Prochaska et al. (2001) 




Q 1224+0037 


1.235 


20.88 


± 


0.06 


>15.11 


Rao, Turnshek & Nestor (2006) 


Meiring et al. (2007) 




Ol 224+0037 


1.267 


20.00 


± 


0.08 


1 4 54+0 09 


Ran Tnrnshpk Rr Np^tnr ("7006^ 


Meiring et al. (2007) 




Q1247+267 


1.223 


19.88 


± 


0.10 


13.97+0.04 


Pettini et al. (1999) 


Pettini etal. (1999) 




Q 1327-206 


0.853 


19.40 


+ 


0.02 


13.90±0.04 


Rao, Turnshek & Nestor (2006) 


Dessauges-Zavadsky et al. 


in prep. 


Ol 331 + 170 


1.776 


21.14 


+ 


0.08 


14 63+0 03 


npccmi crp^-^fivfirl^lrv pt al f7004^ 
J_^CL>c>tlLI^CL> Z-jlL V aUjIVy d til . \ jL\J\J^ J 


Dessauges-Zavadsky et al. 


(2004) 


Ol 151 +11 8 


1.149 


20.23 


± 


0.10 


1 4 74+0 07 


Pettini pt al CI QQQ^ 


Pettini etal. (1999) 




Q1451+123 


2.255 


20.30 


+ 


0.15 


14.33+0.07 


Dessauges-Zavadsky et al (2003) 


Dessauges-Zavadsky et al. 


(2003) 


Q1525+0026 


0.567 


19.78 


± 


0.08 


14.19+0.06 


Rao, Turnshek & Nestor (2006) 


Dessauges-Zavadsky et al. 


in prep. 


Q1622+238 


0.656 


20.40 


± 


0.10 




Rao & Turnshek (2000) 






Q1622+238 


0.891 


19.23 


± 


0.03 




Rao & Turnshek (2000) 






Q1629+120 


0.900 


19.70 


+ 


0.04 




Rao, Turnshek & Nestor (2006) 






Q2048+196 


1.116 


19.26 


+ 


0.08 


>15.22 


Rao, Turnshek & Nestor (2006) 


Dessauges-Zavadsky et al. 


in prep. 


Q2128-123 


0.430 


19.37 


+ 


0.08 


>14.08 


Ledoux, Bergeron & Petitjean (2002) 


Ledoux, Bergeron & Petit. 


ean (2002) 


Q2206-199 


1.920 


20.44 


± 


0.08 


15.30±0.02 


Pettini et al. (2002) 


Pettini et al. (2002) 




Q2230+023 


1.859 


20.00 


± 


0.10 




Dessauges-Zavadsky et al. (2006) 






Q223 1-001 


2.066 


20.53 


+ 


0.08 


14.83±0.03 


Dessauges-Zavadsky et al. (2004) 


Dessauges-Zavadsky et al. 


(2004) 


Q2328+0022 


0.652 


20.32 


± 


0.07 


14.84±0.01 


Rao, Turnshek & Nestor (2006) 


Dessauges-Zavadsky et al. 


in prep. 


Q233 1+0038 


1.141 


20.00 


± 


0.05 


14.38±0.03 


Rao, Turnshek & Nestor (2006) 


Meiring et al. (2007) 




Q2343+125 


2.431 


20.35 


+ 


0.05 


14.66±0.03 


Dessauges-Zavadsky et al. (2004) 


Dessauges-Zavadsky et al. 


(2004) 


Q2348-144 


2.279 


20.59 


+ 


0.08 


13.84±0.05 


Dessauges-Zavadsky et al. (2004) 


Dessauges-Zavadsky et al. 


(2004) 


Q2352-0028 


0.873 


19.18 


± 


0.10 


13.47±0.06 


Rao, Turnshek & Nestor (2006) 


Dessauges-Zavadsky et al. 


in prep. 


Q2352-0028 


1.032 


19.81 


± 


0.14 


14.96±0.04 


Rao, Turnshek & Nestor (2006) 


Dessauges-Zavadsky et al. 


in prep. 


Q2352-0028 


1.246 


19.60 


± 


0.30 


14.28±0.03 


Rao, Turnshek & Nestor (2006) 


Dessauges-Zavadsky et al. 


in prep. 



defining A mm mlix in the Mg n A 2796 A line. We do not appeal to ad- 
ditional lines at this stage, since we are aiming to develop a method 
that can be used with the limited spectral information that is avail- 
able for low redshift absorbers. 



The values of EW and At' that are input into Eqn Q] are gov- 
erned by the choice of /i m i n , max . We have investigated the follow- 
ing possibilities: (i) /l m i n , max defined by 3 a absorption (described 
above and in Ellison 2006); (ii) /t m j n . max defined by the central 90% 
of the Mgn line EW. This is analogous to the technique used on 
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Figure 2. Two examples of fits to high resolution (top panels) and FWHM=5 A (bottom panels) for the sub-DLAs towards Ql 101-264 (left panels) and 
Q2344+0342 (right panels). In all panels the red dotted line shows the damped profile of the best fit to the low resolution data (see red panel labels) and the 
blue dashed line shows the fit by Dessauges-Zavadsky et al. (2003) to the UVES data (blue panel labels). The profiles have been convolved to the resolution 
of the data in each panel. This Figure demonstrates that when the sub-DLA is isolated, the N(H i) determination is excellent, even at moderate resolution. 
However, when hya forest blending becomes significant, there is a tendency to over-estimate the N(H i) from low resolution data. 



unsaturated lines to determine the velocity width of an absorber 
based on optical depth (e.g. Prochaska & Wolfe 1997; Ledoux et 
al. 2006); (iii) /i m i 11>max defined by the central 90% of the velocity 
spread of the Mg n line. As a further experiment, we apply these 
three techniques using the full spread of all Mgn components in 
our calculation of the £>-index, in addition to the original definition 
of Ellison (2006) which uses only the main absorption complex. 
The D-index calculated over the central absorption excludes outly- 
ing components once the continuum has been recovered, whereas 
'full' refers to all absorption components in a given Mg n system. 
In Figure[3]we show how these different approaches yield different 
values for /l nl i n . raax in one of the absorbers in our sample. This Figure 
is reproduced for all of our absorber sample in the online material 
that accompanies this article. We also provide in the online mate- 
rial a table listing all of the EWs and errors for the 3 approaches 
described above for each absorber; Table|4]gives the first four lines 
of the online version as an example. 



In Table [5]wegive the success statistics and optimal D value 
for each techniqu4j. The success rate is defined as the fraction of 
systems with D > A>ptimai that are actually DLAs. The DLA miss 
rate is the fraction ofDLAs with D < D of ,i nwi \. We also tabulate the 
Kolmogorov-Smirnov (KS) probability that the £>-indices of sub- 
DLAs and DLAs for a given technique are drawn from the same 
population. Table [5] shows that all of the techniques yield more or 
less consistent results, with success rates ~ 50 - 55 % and KS prob- 
abilities on the order of a few percent. As an example, we show the 
distribution of D versus N(H i) in Figuref4]when the 90% EW limits 
are applied to the central absorption complex. 



We next consider whether including information from other 
metal lines may improve the use of D to distinguish DLAs and 

3 The D-index statistics which include N(Fe n) are described later in this 
section. 
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Table 4. Equivalent widths and velocity spreads for lambda limits discussed in the text (see online material for the full table). 



QSO 


log N(HI) 


EW 90%EW (A) 


Av 9 o %£ n/(kms ') 


EW 90 %,w(A) 


Av 90 %wi (kms ') 


EW 3ct - (A) 


Av 3[r (km s ') 


Q0823-223 


19.04±0.04 


0.85±0.03 


186.38 


0.93±0.04 


231.30 


0.93±0.04 


206.59 


Q 1206+459 


19.04±0.04 


0.56±0.02 


106.86 


0.63±0.03 


142.48 


0.63±0.03 


115.76 


Q0424-131 


19.04±0.04 


0.25±0.01 


37.53 


0.28±0.02 


65.04 


0.26±0.02 


40.03 


Q0002+051 


19.08±0.04 


0.62±0.03 


95.03 


0.66±0.04 


129.80 


0.63±0.02 


85.76 



Table 5. D-index statistics for /i m j n . max tests 



Absorption range 


A limits 


Optimal D 


DLA success rate (%) 


DLA miss rate (%) 


KS probability (%) 


Central 


3cr 


6.7 


48.6 


19.0 


9.1 


Central 


90% velocity 


5.5 


54.3 


9.5 


0.1 


Central 


90% EW 


7.2 


54.1 


4.8 


0.4 


Full 


3a- 


6.3 


51.5 


19.0 


5.3 


Full 


90% velocity 


5.2 


53.1 


19.0 


2.4 


Full 


90% EW 


7.0 


56.7 


19.0 


1.3 


Central + N(Fen) 


3cr 


6.3 


48.4 


11.8 


20.6 


Central + N(Fen) 


90% velocity 


5.7 


57.7 


11.8 


0.1 


Central + N(Fen) 


90% EW 


7.0 


57.1 


5.9 


0.5 


Full + N(Fen) 


3cr 


5.2 


55.6 


16.7 


1.3 


Full + N(Fen) 


90% velocity 


5.1 


60.0 


16.7 


0.6 


Full + N(Fen) 


90% EW 


5.6 


53.6 


16.7 


1.3 



sub-DLAs. It has been shown (e.g. Meiring et al. 2008; Peroux et 
al. 2008) that sub-DLAs have a tendency to be more metal-rich 
than DLAs. Indeed, there is a lack of high A^H i) high metallic- 
ity systems that empirically manifests itself as an apparent anti- 
correlation between A'(Hi) and [Zn/H] (Boisse et al. 1998; Prant- 
zos & Boissier 2000). Although it has been argued that this is due 
to dust bias, this interpretation is not supported by extinction mea- 
surements of DLAs (Ellison, Hall & Lira 2005). Simply put, the 
observed reddening values derived from DLA samples (Murphy & 
Liske 2004; Ellison et al. 2005; Vladilo, Prochaska & Wolfe 2008) 
are much lower than would be required to impose a 'dust filter'. 
Moreover, it has been argued (Dessauges-Zavadsky et al. 2003) that 
the higher metallicities in sub-DLAs are not a result of ionization 
corrections. As we have shown above, any systematic error in the 
A^H i) of low redshift sub-DLAs determined from low resolution 
spectra is both small, and in the wrong sense to explain the dif- 
ferent metallicity distributions. In the absence of any experimental 
reason (such as dust bias or ionization correction) behind the ele- 
vated metallicities in sub-DLAs, we explore whether information 
on the metallicity of an absorber can be combined with the kine- 
matic description encapsulated in EqnQ]to improve the selection 
success of the D-index. 

In Figure[5]we show the distribution of [Fe/H] versus D-index 
(as defined by the 90% EW range for the central absorption com- 
ponent) for the 49 absorbers for which N(Fen) measurements are 
available in the literature. It can be seen that in the regime of high 
D-index, the DLAs have preferentially lower metallicities than the 
sub-DLAs. For D > 7, the mean [Fe/H] is -0.67 for sub-DLAs and 
-1.43 for DLAs. These mean values are calculated by treating the 
lower limits of the three absorbers with fully saturated Fe n profiles 
as detections. The difference between the mean metallicities could 
therefore be even larger. Although the DLAs in our sample have 
a higher mean redshift than the sub-DLAs (see next section) the 
evolution of DLA metallicity is very mild. Kulkarni et al. (2005) 
show that the DLA metallicity changes by less than 0.3 dex from 
z ~ 0.5 to z ~ 2.5. Of course, we can not use the actual metal- 



licity to fine-tune the definition of the D-index, since this requires 
a measurement of A'(Hi), the very quantity we are hoping to se- 
lect for. We therefore consider the use of N(Fen). In Figure [6] we 
show histograms of log N(Fe n) for the compilation of DLAs and 
sub-DLAs in Dessauges-Zavadsky et al. (in prep.) at z a bs < 1.5 (the 
redshift regime for which we are trying to develop a selection tech- 
nique). The column density of Fen is a quantity that can be mea- 
sured with relative ease from almost any high resolution spectrum 
suitable for D-index determination. There are almost a dozen differ- 
ent Fe ii transitions with a large dynamic range in /-value at wave- 
lengths not far from Mg n. From Figure|6]we can see that the DLAs 
clearly have systematically higher N(Fe n) than the sub-DLAs. The 
mean log N(Fen) (and standard deviation) is 15.15±0.08 for DLAs 
and 14.64±0.07 for sub-DLAs. Although the element zinc is usu- 
ally considered the best indicator of actual metallicity, we empha- 
size here that we are simply using iron from the empirically moti- 
vated viewpoint that its column density is systematically different 
in DLAs and sub-DLAs. 

We factor N(Fe n) into our calculation of the D-index by mul- 
tiplying the value derived from EqnQ]by log N(Fen) and dividing 
by 15 (to re-scale the D-index): 



EW(A) logN(Fell) 
D = . x s - x 1000 



(2) 



Av(km/s) 15 

In Table[5]we give the success and miss rates for this modified 
D-index definition, again calculating its value for the six different 
combinations investigated above and excluding systems which only 
have N(Fen) limits. Figure|7]shows the distribution of D-index ver- 
sus A'(H i) for our sample using the definition of D given in Eqn 
[2] Although the success rates generally improve slightly (by a few 
percent) when N(Fe n) is included, this is not a significant improve- 
ment given the modest-sized sample. The close consistency of suc- 
cess rates with/without N(Fen) included in the calculation is prob- 
ably due to the small relative difference (in the log) between iron 
column densities. It is also worth noting that in Figure[5]there is ap- 
parently no correlation between D-index and [Fe/H] in DLAs, in- 
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Figure 1. Comparison of the best fit iV(Hi) column densities from high 
resolution UVES data and VPFITs to the UVES data after convolution with 
a gaussian of FWHM = 5 A. The column densities are given in Table|5] 




Velocity (km/s) 

Figure 3. The Mg n A 2796 transition for one of the absorbers in our sam- 
ple. The vertical lines show the minimum and maximum velocities used 
to calculate the D-index for the three different possibilities investigated in 
this paper: 90% of the equivalent width (dashed), 90% of the velocity width 
(dot-dashed) and 3 <x absorption significance (dotted). The respective D- 
indices calculated from these 3 techniques are listed in the top right of the 
panel. Mgn spectra are given for the full sample in the online material that 
accompanies this article. 
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Figure 4. The D-index versus A^H i) for the central range of Mg n absorp- 
tion /l m i n ma x as defined by the 90% EW range. The horizontal dashed line 
shows the optimal D-index as given in Tableland the vetrical dashed line 
demarcates DLAs and sub-DLAs. The open points have z- d bs > 1-5 and filled 
points have Zubs < 1.5, see the discussion in Section l4~71 
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Figure 5. The D-index versus [Fe/H] for the central range of Mg n absorp- 
tion /l m i n>raax as defined by the 90% EW range. Solid points denote DLAs, 
open points denote sub-DLAs. The horizontal dashed line represents so- 
lar iron abundance, for which we adopt the value log(N(Fe)/H i) Q = —4.55 
(Asplund et al. 2005). 



dicating that the D-index does not preferentially select high or low 
metallicity DLAs. On the other hand, 4/5 sub-DLAs with D > 8 
have [Fe/H] > -0.5, the highest values in our sample. Therefore, 
using low values of D to select sub-DLAs may miss the highest 
metallicity absorbers. 



4 SUMMARY AND DISCUSSSION 

We have presented a sample of 56 Mgn absorbers with known 
A'(Hi) column densities and investigated techniques to screen for 
DLAs. The D-index combines measures of both Mg n EW and ve- 
locity spread, and we have also investigated incorporating informa- 
tion on the Fe n column density. The results are robust to the various 
parameterizations that we use, with a typical DLA yield of ~ 55%, 
compared to ~ 35% when Mg n EW alone is used (Rao et al. 2006). 
The relative insensitivity to the precise definition of the D-index is 
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z<1.5 
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N(FeII) 

Figure 6. Distribution of Fe n column densities for DLAs (open histogram, 
bold outline) and sub-DLAs (shaded histogram). 




Log N(HI) 

Figure 7. The D-index versus N(H i) for the central range of Mg n absorp- 
tion /l m i 11>max as defined by the 90% EW range. The D-index has been cal- 
culated from Eqn|2] These data points refer to the row in bold in Table \5\ 
The horizontal dashed line shows the optimal D-index as given in Table 
|3]and the vetrical dashed line demarcates DLAs and sub-DLAs. The open 
points have Zubs > 1-5 and filled points have z&s ^ 1-5, see the discussion 
in Section ETTl 

reassuring, in the sense that it demonstrates that its success is un- 
likely to be a fluke of our dataset or choice of parameterization. 

4.1 Using the D-index 

Given the results in Table [5] what is the 'best' definition of D? 
Given the sample size (and hence, rounding errors in the success 
and miss rates), simply selecting the technique with the best statis- 
tics is not necessarily the best choice. Instead, we consider which 
technique might be the most robust against issues such as noise, 
resolution and absorber environment. Any of the methods that use 
the 'full' velocity range will be sensitive to outlying components 
that may have nothing to do with the main absorber, e.g. be associ- 



ated with a satellite galaxy. As discussed in Section l3TI sub-DLAs 
in particular are known to often have companion absorbers. S/N is 
a consideration for techniques that use a 3 cr cut-off", either over the 
full absorption range, or just the central complex. However, Ellison 
(2006) showed that S/N only significantly degrades the efficiency 
of D-index selection in the central absorption component at S/N 
«: 10. The disadvantage of using the 90% velocity range is that 
it is sensitive to small differences (between spectra) of resolution. 
As Ellison (2006) showed, this is also the case for a 3 cr definition. 
Taking these issues into account, and then re-visiting the statistics 
of Table[5] we suggest that the 90% EW over the central absorption 
range (i.e excluding outlying, low EW clouds), should be the most 
robust to these various troubles and yields a high success rate, with 
only one DLA missed. The optimal values of D are very similar for 
the central/90% EW technique whether or not N(Fe n) is included 
- 7.0 and 7.2 respectively. The KS test probability is slightly better 
when Fe n is included. 

One potential bias in the evaluation of the D-index in the cur- 
rent sample is that the redshift distributions of the DLAs and sub- 
DLAs are not the same. The mean DLA redshift is z a t, s = 1.63, 
compared with 1.06 for the sub-DLAs, with very few of the latter 
above z a bs = 1-5. This potential bias can be seen in Figures|4]and 
[7] where the open points show absorbers at z a b s > 1.5. Confining 
our analysis to absorbers only in this low redshift range includes 
only 5 DLAs with Fe n detections, insufficient for robust statistical 
analysis. However, Mshar et al. (2007) have shown that neither the 
EW distribution, nor the velocity spread of strong Mgn systems 
evolves with redshift. These authors do, however, find differences 
in the number and distribution of 'subsystems' between their low 
and high redshift sample. Whilst the low redshift absorbers tend to 
comprise a single strong absorption component with one satellite 
subsystem, there are more subsystems in high redshift absorbers, 
over the same velocity range. The D-index is sensitive to such dif- 
ferences. However, the redshift differences in kinematic structure 
reported by Mshar et al. (2007) and described above can not ex- 
plain the 'success' of the D-index in distinguishing DLAs from 
sub-DLAs. This is because the redshift inequality in our sample 
is such that we have more DLAs at high z, where the tendency 
towards more complex kinematic structure would decrease the D- 
index. However, DLAs are defined by high values of D, i.e. with 
few kinematically extended subsystems. Further evidence against a 
redshift bias in our sample comes from Ledoux et al. (2006) who 
find that, in DLAs, the median velocity width of unsaturated metal 
lines increases with decreasing redshift. This redshift dependence 
again works in the opposite sense to a redshift bias that would ar- 
tificially introduce a D-index dependence in our work. Therefore, 
although it is certainly desirable to repeat the D-index tests for a 
sample of sub-DLAs and DLAs that are matched in redshift, there 
is no empirical evidence to suggest that redshift evolution can ex- 
plain the difference in D. 

In summary, we recommend that the D-index is calculated us- 
ing the 90% EW range over the central absorption complex, incor- 
porating N(Fe it), if available, see Equation^ In Figure|7]we show 
the distribution of N(H i) versus D-index for this choice of parame- 
ters. This Figure is directly comparable to Figure [4] which uses the 
same definition of the D-index, but does not include N(Fe n) in the 
calculation. Based on our recommended definition, the KS proba- 
bility that the D-indices of sub-DLAs and DLAs are drawn from 
the same population is 0.5%. In our sample of 45 absorbers with 
accurate Fen column densities, 57+}^% (16/28) of systems with 
D > 7.0, are DLAs where the error bars are 1 cr values taken from 
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Gehrels (1986). Only 6% (1/17) of the DLAs in the sample of 45 
have D < 7.0. 



4.2 Interpretation of the D-index 

Ellison (2006) suggested that the £>-index may be driven in part by 
the QSO-galaxy impact parameter. For a galaxy composed of nu- 
merous Mgn 'clouds', whose size is of the order of a few kpc (Elli- 
son et al. 2004b), the absorption may appear more patchy when the 
sightline intersects at large impact parameter where the filling fac- 
tor of these clouds is low. A sightline passing close to the galaxy's 
centre will intercept a higher density of absorbing clouds, perhaps 
with a smaller radial velocity distribution. This picture is supported 
by empirical evidence that the EW of Mg n absorbers does corre- 
late with impact parameter (Churchill et al. 2000; Chen & Tinker 
2008) and sub-DLAs do have more distinct absorbing components 
than DLAs (Churchill et al. 2000). 

Another possible source of distinction between the kinematic 
structure of DLAs and sub-DLAs is the presence of galactic out- 
flows. The connection between outflows and a subset of strong 
Mg ii absorbers has been discussed in the literature for several years 
(e.g. Bond et al. 2001; Ellison, Mallen-Ornelas & Sawicki 2003), 
although it remains contentious as a general description (see dis- 
cussions in Bouche et al. 2006 and Chen & Tinker 2008). Nonethe- 
less, galactic outflows are apparently common in high redshift star- 
forming galaxies (e.g. Shapley et al 2003; Weiner et al. 2008) and a 
growing body of empirical evidence is being published which sup- 
ports the connection of strong Mg n absorbers with winds or out- 
flows. For example, Zibetti et al. (2007) find that the mean impact 
parameter of Mg n absorbers in the SDSS is about 50 kpc from a 
L* galaxy. One explanation for this observation is that the absorp- 
tion detected in QSO spectra corresponds to material that is in a 
large extended region, possibly associated with winds. Menard & 
Chelouche (2008) also argue for the connection of Mg n absorbers 
to massive galaxies based on their gas-to-dust ratios, which are an 
order of magnitude higher than values for DLAs (Murphy & Liske 
2004; Ellison et al. 2005; Vladilo et al. 2006). Murphy et al. (2007) 
observed a significant correlation between Mg n EW and metallic- 
ity and concluded that '[there is] a strong link between absorber 
metallicity and the mechanism for producing and dispersing the 
velocity components'. Finally, Kacprzak et al. (2007) have shown 
that there is a correlation between the Mg n EW and the degree of 
asymetry in the host galaxy. However, it does not necessarily flow 
that all strong Mgn absorbers are connected with outflows. Both 
Bouche (2008) and Menard & Chelouche (2008) have recently 
suggested that, when plotted in the ,/V(Hi)-Mgii EW plane, there 
are two separate absorber populations distinguished by metallicity. 
This effect is also seen in Figure [5] where the high D-'mdex ab- 
sorbers are fairly cleanly separated into high metallicity sub-DLAs 
and low metallicity DLAs. Bouche (2008) proposes that the more 
metal-rich, lower A^H i) absorbers at a given Mg n EW are asso- 
ciated with galactic outflows. In turn, this may lead to kinematics 
that are more patchy in velocity space (e.g. Ellison et al. 2003) and 
contribute to the mechanism behind the success of the D-'mdex. 



niche for the £>-index? Individual groups have already used mod- 
erately large datasets for Mg n searches at high resolution. The two 
most recent surveys for strong Mg n absorbers in high resolution 
spectra are Mshar et al. (2007) (56 absorbers with A 2796 EW > 0.3 
A) and Prochter et al. (2006a) (21 absorbers with Mg n A 2796 EW 
> 1 A). Given that the number density of Mg n A 2796 EW > 0.6 A 
absorbers is a factor 2 larger than at 1 A, these two samples would 
already yield over 50 absorbers suitable for D-'mdex screening (not 
accounting for an overlap between the samples). However, these 
samples are small compared to the full archival power that could 
be applied. After a decade of high resolution spectroscopy with 
HIRES on Keck and UVES at the VLT, and a somewhat shorter, 
but extremely productive history with ESI, we estimate that 500 - 
700 QSOs have been observed at moderate-high resolution (e.g. 
Prochaska et al. 2007b). If we assume an average emission red- 
shift z cm =2.5 for 500 QSOs, and maximum wavelength coverage 
of 8500 A, this yields a redshift path of almost 760. The num- 
ber density of strong Mg n absorbers has been robustly determined 
from SDSS spectra (e.g. Nestor et al. 2005; Prochter et al. 2006a). 
For example, Nestor et al. (2005) find that the number density of 
Mgn A 2796 EW > 0.6 A absorbers (a typical EW suitable for 
DLA pre-selection) at z ~ 1 is ~ 0.5. This value is unlikely to 
be affected due to magnitude bias (Ellison et al. 2004a) and can 
therefore also be applied to high resolution spectra. The total num- 
befl of such absorbers that we could therefore expect to find in 
existing moderate-high resolution spectra is therefore approaching 
400. Even though follow-up observations with space telescopes are 
likely to focus on the brightest in this sample, the typical magnitude 
limit for echelle spectroscopy with an 8-m ground-based telescope 
is close to the limit that will be feasible for the Cosmic Origins 
Spectrograph (COS) and the renovated STIS . 

The D-'mdex may also be useful for assessing the nature of ab- 
sorption where follow-up observations of hya are not even possi- 
ble. For example, in the spectra of GRBs. As observatories (and ob- 
servers) fine-tune their rapid follow-up techniques, there is a grow- 
ing database of bright GRBs that have been observed promptly at 
high resolution (e.g. Vreeswijk et al. 2007; Ellison et al. 2007; 
Prochaska et al. 2007a). Mgn absorbers detected in these spec- 
tra that do not have Lyor covered, have no chance of susbsequent 
follow-up due to the rapid fading of the optical afterglow, but their 
likelihood of being DLAs can still be assessed, (e.g. Ellison et al. 
2006; Prochaska et al. 2007a). The application to GRB spectra is 
particularly interesting given the possibility of deep searches for the 
absorbing galaxy after the optical afterglow has faded. Even at low 
redshifts, the study of DLA host galaxies has been challenging with 
QSO impact parameters of 1-2 arcseconds (e.g. Chen & Lanzetta 
2003 and references therein). Identifying probable DLAs in GRB 
spectra opens the door to host galaxy searches to unprecedented 
depths. 

Finally, in light of the recent puzzling result by Prochter et 
al. (2006b) that intervening Mg n absorbers are more numerous to- 
wards GRB sightlines than QSO sightlines, it would also be inter- 
esting to compare the distribution of £>-indices of these two popu- 
lations. 



4.3 Applications of the £>-index 

The £>-index requires moderately high resolution spectra to yield 
accurate pre-selection. Therefore, whilst EW-only DLA pre- 
selection can exploit the truly enormous datasets of, for example, 
the SDSS (e.g. Rao & Turnshek 2000; Rao et al. 2006), what is the 



We note that the sample used for the study in this paper does not have 
much overlap with this archival path length. 
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